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ABSTRACT 

In this Letter we explore the hypothesis that the quasar SDSSJ092712. 65+294344.0 
is hosting a massive black hole binary embedded in a circumbinary disc. The lightest, 
secondary black hole is active, and gas orbiting around it is responsible for the blue- 
shifted broad emission lines with velocity off-set of 2650 km s~^, relative to the galaxy 
rest frame. As the tidal interaction of the binary with the outer disc is expected to 
excavate a gap, the blue-shifted narrow emission lines are consistent with being emitted 
from the low-density inhomogeneous gas of the hollow region. From the observations 
we infer a binary mass ratio q ^ 0.3, a mass for the primary of Mi « 2 x 10^ Mq, 
and a semi-major axis of 0.34 pc, corresponding to an orbital period of 370 years. We 
use the results of cosmological merger trees to estimate the likely-hood of observing 
SDSSJ092712.65-f 294344.0 as recoihng black hole or as a binary. We find that the 
binary hypothesis is preferred being one hundred times more probable than the ejection 
hypothesis. If SDSSJ092712. 65-^294344.0 hosts a binary, it would be the one closest 
massive black hole binary system ever discovered. 

Key words: quasars: individual (SDSSJ092712.65-f 294344.0) — galaxies: kinematics 
and dynamics — galaxies: nuclei — black hole physics 



1 INTRODUCTION 

If massive black holes (MBHs) are ubiquitous in galaxy 
spheroids and galaxies experience multiple mergers dur- 
ing their cosmic assembly, MBH binaries (MBHBs hereon) 
should be common, albeit transient features of galactic 
bulges. Observationally, the paucity of dual active nuclei 
points toward rapid inspiral of the MBHs down to parsec 
scales where they form a gravitationally bound pair (Mayer 
et al. 2007). After birth, the MBHB further hardens under 
the action of gas/dynamical and/or stellar torques (Begel- 
man, Blandford & Rees 1980) and depending on their effi- 
ciency (Milosavljevic & Merritt 2001; Armitage & Natarajan 
2002; Berczik et al. 2006; Sesana et al. 2007; Perets et al. 
2007; Colpi et al. 2007) the binary may remain in its hard 
state for a relatively long time, before coalescing under the 
action of gravitational waves. 

MBHBs can coalesce shortly, and a clear prediction of 
general relativity is that the relic MBH receives a kick in 
response to the asymmetric pattern of gravitational waves 
emitted just prior coalescence. For slowly spinning MBHs 
recoil velocities of a few x 100 km s~^ were predicted, but re- 
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cent numerical advances in general relativity have indicated 
that velocities as extreme as '--^ 4000 km s~^ can be attained 
for maximally rotating Kerr MBHs with spin vectors lying 
in the orbital plane (Schnittman & Buonanno 2007; Lousto 
& Zlochower and references therein) . These velocities far ex- 
ceed the escape speed even from the brightest galaxies and 
for this reason MBHs receding from their parent host have 
been considered as electromagnetic targets of a coalescence 
event (Devecchi et al. 2008; Loeb 2007; Volonteri & Madau 
2008). 

Recently, the quasar SDSS J092712. 65-1-294344.0 (here- 
after SDSSJ0927; Adelman-McCarthy et al. 2007) has been 
considered as first candidate of a recoiling MBH by Ko- 
mossa, Zhou & Lu (2008, KZL). SDSSJ0927 exhibits, in its 
optical spectrum, two distinct sets of lines offset by 2650 
km s~^ relative to each other. The first set of narrow emis- 
sion lines (NELs) has redshift z = 0.713 (identified by KZL 
as the redshift of the galaxy that hosted the MBH at its 
birth). We will refer to this set that is typical of an AGN as 
rest-frame NELs (rf-NELs hereon; the red NELS in KZL). 
The rf-NELs do not show any evidence of ionization strat- 
ification (as in few other quasars), possibly because these 
lines are not resolved in the SDSS spectrum. The second set 
of lines comprises two blue-shifted systems featuring dif- 
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ferent FWHM: the broad Mg II and Balmer emission lines 
(b-BELs hereon) with FWHM « 4000 km s"\ and the high- 
ionization narrow lines (b-NELs) with FWHM « 460-2000 
km s~^, both consistent with a redshift z — 0.698. The b- 
BELs and b-NELs refer to a line system in coherent motion 
relative to the rest-frame of the host, with a light-of-sight ve- 
locity of 2650 km s~^. In KZL, the permitted b-BELs result 
from a broad line emission system gravitationally bound to 
the recoiling MBH, while b-NELs remains problematic for 
the reasons explained in Section [2] 

The aim of this Letter is to suggest an alternative in- 
terpretation to SDSSJ0927: that SDSSJ0927 is an unequal- 
mass MBHB embedded in a circumbinary gaseous disc at the 
centre of its host galaxy. In Section |3] and |4] we describe 
the main properties of the model that are requested to re- 
produce the observed properties of b-BELs. In Section Owe 
discuss the possible origin of b-NELs and of their Doppler 
shift. In Section |B] we present our conclusions. 



2 LIMITATIONS OF THE EJECTION 
HYPOTHESIS 

The origin of the b-NELs and the likelihood of observing a 
recoiling MBH with a speed close to its maximum are the 
main concern of this Section. 

First, we notice that the recoil hypothesis by KZL seems 
unable to explain all properties of the b-NELs observed in 
SDSSJ0927. The FWHM of the b-NELs is larger than the 
FWHM of the rf-NELs and of typical FWHMs of NELs 
observed in AGNs. By contrast, the FWHM of b-NELs is 
not large enough to assure that the emitting gas is bound 
to the recoiling MBH, although this should be the case as 
the b-NELs have the same Doppler shift as the b-BELs. 
KZL suggest that the b-NELs are emitted by gas shocked 
in the interaction of the MBH with the interstellar medium. 
In this case however, there is no reason why the speed of the 
shocked gas should equal that of the MBH (Devecchi et al. 
2008). Alternatively, KZL consider that b-NELs are emitted 
in the outer part of the accretion disc with substantial radial 
spreading and/or in an unspecified kind of outflow. This 
last possibility was discussed in details by Bogdanovic et al. 
(2008). They find that mass outflows strongly reduce the 
density of the inner regions of the accretion disc implying 
a low rate of gas accretion onto the MBH. In this case the 
luminosity would be much weaker than observed, unless the 
kicked MBH is more massive than 5 x 10^ Mq. 

Second, we remark that the probability of observing 
an ejected MBH with the properties of SDSSJ0927 depends 
sensitively on the specific frequency associated to MBHB 
mergers able to form a remnant with the required large mass 
(Mbh), and on the probability for the remnant to receive 
such a high recoil velocity, i.e. v > 2500 km s~^ (we conser- 
vatively adopt this lower limit to the recoil velocity, which 
provides an upper limit to the detection probability). The 
latter constraint strongly depends on the MBH spin mag- 
nitudes and orientations with respect to the orbital plane, 
and on the mass ratio of the binary, q = M2/M1 ^ 1 as il- 
lustrated in Fig. 1. We model the recoil using the two avail- 
able fits to the recoil velocity based on the latest numerical 
relativity results (Lousto & Zlochower 2009). Non-spinning 
holes always recoil with velocities below 200kms~^, and the 
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Figure 1. Bottom panel; recoil velocity as a function of the bi- 
nary mass ratio q = M2/M1 , for maximally Kerr holes with a = 1. 
Long-dashed curve: maximum recoil. Short-dashed curve; min- 
imum recoil. Solid curve: mean recoil, with 1— ct scatter. In all 
cases the recoil velocity is averaged over 10^ orbital configura- 
tions for each q. Top panel: probability for recoils faster than 
V > 2500 km s~^ as a function of q, for three representative spin 
values (both holes are assumed to have the same spin) , assuming 
isotropic orbital parameters. No recoil with v > 2500 km s~^ is 
possible for MBH spin parameters d < 0.65. 



recoil increases with increasing spin. In general the recoil is 
maximized for spin vectors lying in the orbital plane. 

We can estimate the probability of detection of an ejec- 
tion with V > 2500 km s~^ by convolving the probability 
of large recoils (shown in Fig. 1) with the merger rate of 
MBHBs. We can compute the number of MBH coalescences 
observable across the entire sky, based on the merger rates 
by Volonteri, Haardt & Madau (2003; Volonteri & Madau 
2008). For a statistical sample of MBH mergers, our models 
provide the masses of the merging MBHs and the redshift 
of the event. We select coalescences in the redshift range 
0.5 < 2 < 1 (compatible with the redshift of SDSSJ0927), 
further imposing a total binary mass Mbh > 5 x 10* M© (i.e. 
the mass of the hypothetical recoiling MBH as discussed in 
KZL2008). Marginalising over MBH masses and q, the high- 
est rate of successful recoils, v > 2500 km s^^, is found by 
imposing maximal values for the dimensionless spin parame- 
ters, a = 1 for both holes; the spins initially lie in the orbital 
plane pointing in opposite directions. The rate is 5 x 10~® 
yr~^. The requirement that the b-NELs come from a spread- 
ing accretion disc however implies Mbh > 5 x 10^ Mq . For 
these masses, the rate of events with v > 2500 km s~^ de- 
creases by one order of magnitude. This set of assumptions, 
maximizing the recoil over the orbital parameters, is mostly 
dependent on the mass ratio q. We need mass ratios 5 ~ 1 
to obtain recoils above 10^ kms~^. Equal mass mergers are 
however rare (e.g., Volonteri & Madau 2008; Sesana et al. 
2005); the distribution of mass ratios for Mbh ^ 5 x 10® 
M0 at 0.5 < z < 1 peaks at q = 0.3 and decreases steeply at 
larger q. Assuming MBHs with equal spins, a = 1, isotropi- 
cally distributed, and Mbh > 5 x 10* Mq, the expected rate 
decreases to less than 10~^ yr~^- The rate of ejections de- 
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the low density of the gas. Since M2 and the gas orbiting in 
the gap at comparable annuli are subject to the same gravi- 
tational potential of Mi , both move with approximately the 
same Keplerian velocity. Accordingly, the same blue-shift 
for the b-BELs and b-NELs can be explained if M2 emits 
an anisotropic ionizing radiation normal to the plane of the 
discs so that the ionizing photons interact only with the gas 
in the gap near M2. This is naturally achieved if, e.g., the 
accretion flow around M2 is a standard Shakura-Sunyaev 
accretion disc. Indeed, cosine-like behaviour of the thermal 
ionizing flux, coupled to the small solid angle subtended by 
the disc assures that the ionization parameter of the gap gas 
not in the immediate vicinity of M2 is relatively low. 



Figure 2. Sketch of the MBHB circumbinary disc structure in 
the nucleus of SDSSJ0927. 



creases drastically with decreasing spin magnitude, and it is 
zero for any orbital-spin conflguration if a < 0.65. 



3 THE MBH BINARY MODEL 

In our model, SDSSJ0927 hosts a MBHB surrounded by a 
circumbinary thin disc feeding the secondary MBH. In the 
tidal interaction of the binary with the gas a gap, i.e., a low 
density region surrounding the binary, opens. 

It is known that, for MBHBs with q ~ 1, A/2 is massive 
enough to perturb strongly the surrounding disc, creating 
over few orbital periods an annular low density region, com- 
monly referred as "gap" . During this phase most of the inner 
disc is accreted by the primary (Syer & Clarke 1995; Ivanov 
et al. 1999, Dotti et al. 2006b), preventing long lasting ac- 
cretion events onto Mi. 

As shown in Section 4, in order to explain the spectrum 
of SDSSJ0927, q ~ 0.3. Given this constraint, we expect that 
the tidal interaction of the binary with the viscous disc will 
leave Mi in a state of low activity. M2, remaining in closer 
gravitational contact with the inner rim of the outer disc, 
moves at low speed relative to the higher density reservoir of 
gas, thus sweeping up an accretion stream that inflows from 
the edge of the gap. This continuous refueling supplies an 
accretion disc around AI2, and allows for luminous accretion 
events during the orbital decay of M2 (Hayasaki, Mineshige, 
Sudou 2007; Cuadra et al. 2008). The position of the inner 
edge of the outer circumbinary disc for q ~ 0.3. (of interest 
for the scope of the paper) is known to lie between ~ 1.5 — 2a, 
where a is the semi-major axis of the MBHB. As for large q 
migration occurs on timescales much longer than the viscous 
time in the outer disc (Ivanov et al. 1999), this binary can 
be long-lived and we work under this hypothesis. A sketch 
of the model is shown in Fig. [21 

We will consider M2 the only accreting MBH of the bi- 
nary. In this model, the rf-NELs are emitted from the "stan- 
dard" large-scale narrow line region that extends around the 
binary for ~ 1 kpc, while the b-BELs are produced in the 
broad line region gravitationally bound to M2, and can be 
blue-shifted or red-shifted depending on the orbital phase 
of the secondary. The b-NELs come from a region of the 
gap near M2, where forbidden lines are emitted because of 



4 ORBITAL PARAMETERS 

In this section we use some of the observed properties of 
SDSSJ0927 (the blue-shift and the FWHM of the b-BEL, 
and the monochromatic luminosity at 5100 A) to constrain 
the dynamical properties of the hypothetical MBHB. The 
compatibility of the model with the other observed spectral 
features (in particular the blue-shift and the flux of the b- 
NELs) is discussed in Section[5] Given the large blue-shift of 
the observed b-NELs and b-BELs, corresponding to a high 
velocity of the accreting MBH, we assume that the active 
MBH is M2, the lightest hole, and that its mass is 6 x 10* 
M0 as inferred in KZL2008. We will further assume that the 
binary is moving on a (quasi) circular orbit. 

The flrst constraint to be verifled is that the observed 
broad line region responsible for the b-BELs is bound to M2, 
and is not tidally stripped by Mi . This condition can be re- 
written as i?BLR < R2, where -Rblr ~ 0.1 pc (KZL2008) is 
the radius of the BLR and R2 is the Roche lobe radius of 
Al2. We estimate R2 using the approximation in Eggleton 
(1983), and impose: 



Rblr, ^ R2 
a a 
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0.6g2/3 +ln(H-gi/3) 

where a is the semi-major axis of the MBHB. 

The observed velocity, Uobs, of Af2 depends on Mi, a, 
and the orientation of the plane of the binary in the sky: 



■Uobs = 1>2 Sm I COS <p — 



GM2 

aq{l + q) 



(2) 



where V2 is the orbital velocity of A/2 relative to the centre of 
mass, i G [0°,90°] is the angle between the line of sight and 
the direction normal to the orbital plane, and £ [0°, 360°] 
is the phase of the orbit. From eqs. [T] and [2l and assuming 
cos (f) = 1 in order to minimize A/i , we obtain 



GM2 



. 2 • . 1/3 

sm I > q ' 



+ 
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0.49 



l/3x 



(3) 



The fraction on the left hand side of eq.[3]is constrained 
by the observations while the right hand side is a mono- 
tonically increasing function of q. Given the inclination i, 
eq. [3] defines a maximum q, and a corresponding minimum 
Ml. We focus on the minimum mass of Mi, obtained for 
the largest value of sini, that is fltting all the observational 
constraints. Assuming that the unification model for AGNs 
(Urry & Padovani 1995) remains valid for this peculiar ob- 
ject, the observation of the b-BELs limits the maximum 
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value of the inclination of the orbital plane to imax ~ 40°, 
which implies q « 0.35, a ^ 0.34 pc, and Mi ^ 1.7 x lO^M©, 
value consistent with the observed high mass end of the 
MBH mass function. These values correspond to an or- 
bital period of 370 yr. The time-scale for two MBHs in 
SDSSJ0927 to reach final coalescence due to gravitational 
wave emission is 3 x 10^ yr, shorter than the Hubble time. 
If the binary model will be confirmed, this is proof of the 
existence of a dynamical process able to overcome the so 
called "last parsec problem" and drive MBHBs to the final 
coalescence. 

Using the FWHM of [OIII] as a tracer of the stellar ve- 
locity dispersion, a,, and assuming the Mbh — c* relation 
(Ferrarese & Ford, 2005 and references therein), we obtain 
for Ml a mass ~ 3 orders of magnitude smaller than our es- 
timate. This would be in contrast with the high luminosity 
of SDSSJ0927, which would imply extremely super Edding- 
ton luminosity (~ lO'^ Eddington). On the other hand, in 
the MBHB scenario, the dynamics of gas may be not repre- 
sentative of the galactic potential, due to the recent galaxy 
merger. 

To estimate the likelihood of the detection of a MBHB 
with the properties inferred by our model, we use the 
same theoretical merger-tree for cosmological MBHB used 
in Section [21 We again impose limits on the redshift range 
(0.5 < 2 < 1), but now the binary model constrains the 
masses of the binary members, rather than the mass of the 
remnant. We therefore compute the probability that an ob- 
server at z — detects binaries with Mi > 1.5 x 10® M0 
and M2 > 5 X 10** M0 across the entire sky. The observable 
merger rate corresponds to ~ 10~* binaries per year, a value 
20 times larger than the expected rate of ejections obtained 
using the most favourable conditions for high kick velocity, 
as discussed in Section (2] 



5 NARROW EMISSION FROM THE GAP 

Our model suggests that the b-NELs are produced from the 
low-density gas in the gap region near M2 . This assumption 
explains the blue-shift of these lines, as we discussed in Sec- 
tion [3] For consistency, we need to constrain the properties 
of the gas in the gap, and check if its physical status is con- 
sistent with the production of b-NELs and with the energy 
flux observed at these frequencies. 

We assume that the gas in the gap is mostly ionized, so 
that we can derive the gas properties from the luminosity of 
the Balmer lines. Under the hypothesis of isotropic emission 
the total luminosity of the H/3 line is 

Lh0 = ^tvDIFup = huni3npneau0Vgs,p, (4) 

where Dl ~ 4127 Mpc is the luminosity distance of 
SDSSJ0927, Fhp = 1.9 x 10"" erg s"^ cm'^ is the flux 
of the b-NEL component of the H/3 line (that we obtained 
directly by fitting the rest-frame, de-reddened SDSS spec- 
trum), h is the Planck constant, uup is the frequency of the 
H/3 photons, Up and are the densities of protons and elec- 
trons respectively, and an/s is the recombination coefficient. 
In order to estimate we need to constrain the volume of 
the emitting region, Vgap, that we assume to be half of the 
total volume of the gap in which M2 resides (see discussion 
in Section [3]) 



Fgap = 27r / r^^dr, (5) 
Jo ^ 

where H{r) is the thickness of the low density gaseous struc- 
ture at a radius r and -Rgap ~ 2a. We use the standard 
assumption that H{r)/r = Cs/«Kepi, where Cs is the local 
sound speed of the plasma and VKepi is the Keplerian ve- 
locity of the gas due to the potential well of Mi. In this 
first estimate we assume T ~ 10* K, that corresponds to 
H{r)/r ~ 5 X 10^'^. Replacing this ratio in eq. [5]and using 
eq. Uwe obtain ne « 8 x 10^ cm^^. 

A difi'erent estimate for rie can be obtained assuming 
that the vertical support of the low density gas in the gap is 
due to a supersonic turbulent motion of the plasma of the 
order of 1/2 the average FWHM of the b-NELs. In this sec- 
ond estimate we assume Uturb ~ 200 km s~^, corresponding 
to ne ~ 2 X 10^ cm"''. We emphasize that these estimates 
must be considered as a lower limit, because the gas can 
be clumpy (so that only a fraction of the irradiated gas will 
contribute consistently to the H/3 emission), and because the 
anisotropic radiation produced in the accretion disc of M2 
can irradiate less than 1/2 of the gap. 

Both estimates of rie are much smaller than the den- 
sity expected for a standard a-disc around Mi at r ~ a, 
indicating the presence of a gap, confirming our expecta- 
tions. An additional, independent constraint on rie can be 
obtained from the flux ratios of different Oxygen lines, in 
particular the [OII]3727 doublet and the [OIII]4363, 4959, 5007 
lines (see e.g. Osterbrock & Ferland 2006 and references 
therein). The ratio between the two components of the [OH] 
doublet cannot be employed here, since the lines are practi- 
cally unresolved. We restrict our analysis to the [OIII] lines. 
Since the b-[OIII]5007 and the rf-[OIII]4959 are blended, we 
assume a factor of 3 between the fluxes of the [OIII] 5007 
and [OIII]4959 lines in both systems (Osterbrock & Ferland 
2006). We compare the observed line flux ratios to the ex- 
pectations for a gas in photo-ionization equilibrium, which 
yields Te = 10* - 2 x 10* K. For Te = 10* K, we obtain 
Ue — (2.5^0 5) X 10*" cm~^. This density is fully consistent 
with the predictions from the H/3 luminosity. Increasing Te, 
Ue decreases. At the maximum temperature allowed by the 
photo-ionization equilibrium, 2 x 10* K, the b-system has 
Tie = (4ti) X 10^ cm"^. Typically [OIII] lines are produced 
in regions with much lower density (narrow line regions, 
rie ~ 10^ cm"^). SDSSJ0927 is indeed confirmed to be a 
peculiar system. 

We check our results by calculating if the production 
of a forbidden line such as [OIII] in such high density re- 
gions is self-consistent. We therefore consider the contribu- 
tion of both collisional excitations and de-excitations when 
computing the luminosity of [OIII] lines. This calculation 
depends on the volume of the emitting region, on rie, and 
on the density of [OIII] ions, n[oin]- Via detailed equilib- 
rium balance, with the temperatures and electron densities 
derived above, we find n[oin] — 2-4 orders of magnitude 
lower than Tie. This ratio is fully consistent with standard 
metallicities in quasars hosts. 

This consideration can explain the presence of both H/3 
and forbidden lines, but the [OH] doublet. [OH] emission 
is not consistent with such high densities, because the [OH] 
critical density for collisional de-excitation is Ue < 10* cm"''. 
The detection of [OH] lines suggests that the gap is filled 
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with an inhomogeneous, cloudy medium. The source of this 
inhomogeneity may be the competition between the gas in- 
flow from the circumbinary disc and the periodic perturba- 
tion exerted by the MBHB (see, e.g. Hayasaki, Mineshige, 
Sudou 2007; Cuadra et al. 2008; Bogdanovic et al. 2008). In 
an inhomogeneously filled gap, H/3 lines are emitted in the 
highest density regions (Tie > 2 X 10^ cm-3), while [OIII] 
and [on] lines are produced by intermediate (3 x 10^ cm~'^ 
< rie < 3.5 X 10® cm"^) and low (n^ < 10'' cm"^) density 
regions, respectively. 



6 DISCUSSION AND CONCLUSIONS 

The binary hypothesis for SDSSJ0927 relies on a simple set 
of assumptions. We can relax the assumption of circular or- 
bits, and assume a non zero eccentricity (e < 0.5), expected 
to be present after gap opening (Armitage & Natarajan 
2005; Dotti et al. 2006a, 2007; Cuadra et al. 2008). Eccentric- 
ities of this magnitude are still compatible with the binary 
model. Larger values are instead unlikely in SDSSJ0927: if 
the binary has a large eccentricity and M2 is near apocentre, 
an extremely massive Mi is needed to justify such a large 
orbital velocity. On the other hand. Shields et al. (2009) set 
an upper limit on the change of the Doppler shift of the b- 
BEL system (< 10 km/s/year). In the binary scenario, this 
prevents M2 to be caught too close to the pericentre. We 
note that assuming circular orbit and the set of MBHB pa- 
rameters discussed above, the b-BEL system moves only of 
4.2 km/s over three years, corresponding to an average shift 
of « 1.4 km/s/year. 

SDSSJ0927 is clearly exceptional. Either it is the first 
example of a MBH ejected from its nucleus due to the grav- 
itational recoil, or it harbours one of the smallest separation 
MBHBs ever detected. Only the possible binary in OJ287 
could have closer MBHs (~ 4.5 x 10"^ pc, Valtonen 2007). 
SDSSJ0927 would provide clear evidence that MBHBs can 
indeed reach coalescence in less than a Hubble time, and that 
Nature is able to solve the last-parsec problem. If the model 
discussed in this letter is confirmed, SDSSJ0927 will repre- 
sent the first detection ever of emission from a circumbinary 
gap around a MBHB. 

The proof of existence of binaries with such a small 
orbital separation is of extreme importance for gravitational 
wave experiments such as LISA (Bender et al. 1994) , and the 
Pulsar Timing Arrays (PTA, Sesana et al. 2008). Since our 
theoretical expectations for a binary with the characteristics 
of SDSSJ0927 are of order 10~* (binary hypothesis) down 
to 10^® (ejection hypothesis) per year, either we have been 
extraordinarily lucky, or MBHBs are more common than 
predicted, opening exciting possibilities for LISA and PTA. 
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